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Abstract 

This paper is an exposition of the relationship between Witten's functional 
integral and Vassiliev invariants. 

1 Introduction 

This paper shows how the Kontsevich Integrals, giving Vassiliev invariants in 
knot theory, arise naturally in the perturbative expansion of Witten's func- 
tional integral. The paper is a sequel to |26fl . Since the writing of |2(J I 



became aware of the work of Labastida and Perez [|37| on this same subject. 
Their work comes to an identical conclusion, interpreting the Kontsevich in- 
tegrals in terms of the light-cone gauge and thereby extending the original 
work of Frohlich and King |nj . The purpose of this paper is to give an expo- 
sition of these relationships and to introduce diagrammatic techniques that 
illuminate the connections. In particular, we use a diagrammatic operator 
method that is useful both for Vassiliev invariants and for relations of this 
subject with the quantum gravity formalism of Ashtekar, Smolin and Rovelli 
[ fy]] . An aspect that this paper does not treat is the perturbation expansion 
via three-space integrals leading to Vassiliev invariants as in 0. See also [§]. 
Nor do we deal with the combinatorial reformulation of Vassiliev invariants 
that proceeds from the Kontsevich integrals as in ||. 



The paper is divided into three sections. Section 2 discusses Vassiliev 
invariants and invariants of rigid vertex graphs. The section three on the 
functional integral introduces the basic formalism and shows how the func- 
tional integral is related directly to Vassiliev invariants. In this section we 
also show how our formalism works for the loop transform of Ashtekar,Smolin 
and Rovelli. Finally section 4 shows how the Kontsevich integral arises in the 
perturbative expansion of Witten's integral in the axial gauge. One feature 
of section 4 is a new and simplified calculation of the necessary correlation 
functions by using the complex numbers and the two-dimensional Laplacian. 
We show how the Kontsevich integrals are the Feynman integrals for this 
theory. 
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2 Vassiliev Invariants and Invariants of Rigid 
Vertex Graphs 

If V(K) is a (Laurent polynomial valued, or more generally - commutative 
ring valued) invariant of knots, then it can be naturally extended to an 



invariant of rigid vertex graphs |2l] by defining the invariant of graphs in 
terms of the knot invariant via an unfolding of the vertex. That is, we can 
regard the vertex as a "black box" and replace it by any tangle of our choice. 
Rigid vertex motions of the graph preserve the contents of the black box, 
and hence implicate ambient isotopies of the link obtained by replacing the 
black box by its contents. Invariants of knots and links that are evaluated 
on these replacements are then automatically rigid vertex invariants of the 
corresponding graphs. If we set up a collection of multiple replacements at 
the vertices with standard conventions for the insertions of the tangles, then 
a summation over all possible replacements can lead to a graph invariant 
with new coefficients corresponding to the different replacements. In this 
way each invariant of knots and links implicates a large collection of graph 
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invariants. See 



The simplest tangle replacements for a 4-valent vertex are the two cross- 
ings, positive and negative, and the oriented smoothing. Let V(K) be any 
invariant of knots and links. Extend V to the category of rigid vertex em- 
beddings of 4-valent graphs by the formula 

V(K m ) = aV{K + ) + W(K_) + cV(K ) 

where K + denotes a knot diagram K with a specific choice of positive cross- 
ing, K_ denotes a diagram identical to the first with the positive crossing 
replaced by a negative crossing and denotes a diagram identical to the 
first with the positive crossing replaced by a graphical node. 

This formula means that we define V(G) for an embedded 4-valent graph 
G by taking the sum 

V(G) =J2a i+{S) b l - {S) c io{s) V(S) 
s 

with the summation over all knots and links S obtained from G by replacing 
a node of G with either a crossing of positive or negative type, or with a 
smoothing of the crossing that replaces it by a planar embedding of non- 
touching segments (denoted 0). It is not hard to see that if V(K) is an 
ambient isotopy invariant of knots, then, this extension is an rigid vertex 
isotopy invariant of graphs. In rigid vertex isotopy the cyclic order at the 
vertex is preserved, so that the vertex behaves like a rigid disk with flexible 
strings attached to it at specific points. 

There is a rich class of graph invariants that can be studied in this manner. 
The Vassiliev Invariants |46|] , [|7|] , constitute the important special case of 
these graph invariants where a = +1, b = —1 and c = 0. Thus V(G) is a 
Vassiliev invariant if 

V(K t ) = V(K + )-V(K_). 

Call this formula the exchange identity for the Vassiliev invariant V. See 
Figure 1 

V is said to be of finite type k if V(G) = whenever \G\ > k where \G\ 
denotes the number of (4-valent) nodes in the graph G. The notion of finite 
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(K|*) (K| + ) (KI-) 

v(K|*) = v(K|+) - v(Kl-) 



Figure 1 — Exchange Identity for Vassiliev Invariants 
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type is of extraordinary significance in studying these invariants. One reason 
for this is the following basic Lemma. 



Lemma. If a graph G has exactly k nodes, then the value of a Vassiliev 
invariant V}~ of type k on G, v k {G), is independent of the embedding of G. 

Proof. The different embeddings of G can be represented by link diagrams 
with some of the 4-valent vertices in the diagram corresponding to the nodes 
of G. It suffices to show that the value of Vk{G) is unchanged under switching 
of a crossing. However, the exchange identity for v k shows that this difference 
is equal to the evaluation of v k on a graph with k + 1 nodes and hence is 
equal to zero. This completes the proof.// 

The upshot of this Lemma is that Vassiliev invariants of type k are in- 
timately involved with certain abstract evaluations of graphs with k nodes. 
In fact, there are restrictions (the four-term relations) on these evaluations 
demanded by the topology and it follows from results of Kontsevich that 
such abstract evaluations actually determine the invariants. The knot invari- 
ants derived from classical Lie algebras are all built from Vassiliev invariants 



of finite type. All this is directly related to Witten's functional integral 49 



In the next few figures we illustrate some of these main points. In Figure 
2 we show how one associates a so-called chord diagram to represent the 
abstract graph associated with an embedded graph. The chord diagram is a 
circle with arcs connecting those points on the circle that are welded to form 
the corresponding graph. In Figure 3 we illustrate how the four-term relation 
is a consequence of topological invariance. In Figure 4 we show how the four 
term relation is a consequence of the abstract pattern of the commutator 
identity for a matrix Lie algebra. This shows that the four term relation 
is directly related to a categorical generalisation of Lie algebras. Figure 5 
illustrates how the weights are assigned to the chord diagrams in the Lie 
algebra case - by inserting Lie algebra matrices into the circle and taking a 
trace of a sum of matrix products. 
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Figure 3 — The Four Term Relation from Topology 
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Figure 4 — The Four Term Relation from Categorical Lie Algebra 
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Figure 5 — Calculating Lie Algebra Weights 



Vassiliev Invariants and Witten's Functional 
Integral 



In [49| Edward Witten proposed a formulation of a class of 3-manifold in- 
variants as generalized Feynman integrals taking the form Z(M) where 



Z(M) 



DAe (ik/4n)S(M,A)_ 



Here M denotes a 3-manifold without boundary and A is a gauge field (also 
called a gauge potential or gauge connection) defined on M. The gauge field 
is a one-form on a trivial G-bundle over M with values in a representation of 
the Lie algebra of G. The group G corresponding to this Lie algebra is said 
to be the gauge group. In this integral the action S(M, A) is taken to be 
the integral over M of the trace of the Chern-Simons three-form A A dA + 
(2/3)^4 A A A A. (The product is the wedge product of differential forms.) 

Z[M) integrates over all gauge fields modulo gauge equivalence (See || 
for a discussion of the definition and meaning of gauge equivalence.) 

The formalism and internal logic of Witten's integral supports the exis- 
tence of a large class of topological invariants of 3-manifolds and associated 
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invariants of knots and links in these manifolds. 



The invariants associated with this integral have been given rigorous com- 
binatorial descriptions |38| , [j|5| , |28j , , [fiTthPH, but questions and conjec- 
tures arising from the integral formulation are still outstanding. (See for 
example H, , |13|| , [p^],p9|, flOfl .) Specific conjectures about this integral 



take the form of just how it implicates invariants of links and 3-manifolds, 
and how these invariants behave in certain limits of the coupling constant 
k in the integral. Many conjectures of this sort can be verified through the 
combinatorial models. On the other hand, the really outstanding conjecture 
about the integral is that it exists! At the present time there is no measure 
theory or generalization of measure theory that supports it. Here is a formal 
structure of great beauty. It is also a structure whose consequences can be 
verified by a remarkable variety of alternative means. 

We now look at the formalism of the Witten integral in more detail and 
see how it implicates invariants of knots and links corresponding to each 
classical Lie algebra. In order to accomplish this task, we need to introduce 
the Wilson loop. The Wilson loop is an exponentiated version of integrating 
the gauge field along a loop K in three space that we take to be an embedding 
(knot) or a curve with transversal self-intersections. For this discussion, the 
Wilson loop will be denoted by the notation Wk(A) =< K\A > to denote 
the dependence on the loop K and the field A. It is usually indicated by the 
symbolism tr(Pe^K A ) . Thus 

W K (A) =< K\A >= tr(Pe$K A ). 

Here the P denotes path ordered integration - we are integrating and expo- 
nentiating matrix valued functions, and so must keep track of the order of 
the operations. The symbol tr denotes the trace of the resulting matrix. 

With the help of the Wilson loop functional on knots and links, Witten 
writes down a functional integral for link invariants in a 3-manifold M: 



Z(M, K)= J DAe w ^ )s(M > A Hr(Pe$K 
= J DAe {ik/4w)s < K\A > . 
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Here S(M, A) is the Chern-Simons Lagrangian, as in the previous discussion. 
We abbreviate S(M, A) as S and write < K\A > for the Wilson loop. Unless 
otherwise mentioned, the manifold M will be the three-dimensional sphere 
S 3 

An analysis of the formalism of this functional integral reveals quite a bit 
about its role in knot theory. This analysis depends upon key facts relating 
the curvature of the gauge field to both the Wilson loop and the Chern- 
Simons Lagrangian. The idea for using the curvature in this way is due to 
Lee Smolin (See also To this end, let us recall the local coordinate 

structure of the gauge field A(x), where a; is a point in three-space. We can 
write A(x) = A^.(x)T a dx k where the index a ranges from 1 to m with the 
Lie algebra basis {T 1; T 2 , T 3 , T m }. The index k goes from 1 to 3. For each 
choice of a and k, A^.(x) is a smooth function defined on three-space. In 
A{x) we sum over the values of repeated indices. The Lie algebra generators 
T a are matrices corresponding to a given representation of the Lie algebra of 
the gauge group G. We assume some properties of these matrices as follows: 

1. [T a ,7V] = if abc T c where [x,y] = xy — yx , and f abc (the matrix of struc- 
ture constants) is totally antisymmetric. There is summation over repeated 
indices. 

2. tr(T a Tb) = 8 a b/2 where 5 a b is the Kronecker delta (5 a b = 1 if a = 6 and 
zero otherwise). 

We also assume some facts about curvature. (The reader may enjoy 
comparing with the exposition in |23| . But note the difference of conventions 
on the use of i in the Wilson loops and curvature definitions.) The first fact 
is the relation of Wilson loops and curvature for small loops: 

Fact 1. The result of evaluating a Wilson loop about a very small planar 
circle around a point x is proportional to the area enclosed by this circle times 
the corresponding value of the curvature tensor of the gauge field evaluated 
at x. The curvature tensor is written 

F r a s (x)T a dx r dy s . 

It is the local coordinate expression of F — dA + A A A. 
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Application of Fact 1. Consider a given Wilson line < K\S >. Ask how 
its value will change if it is deformed infinitesimally in the neighborhood of a 
point x on the line. Approximate the change according to Fact 1, and regard 
the point x as the place of curvature evaluation. Let 5 < K\A > denote the 
change in the value of the line. 5 < K\A > is given by the formula 

5 < K\A >= dx r dx s F r a s (x)T a < K\A > . 

This is the first order approximation to the change in the Wilson line. 

In this formula it is understood that the Lie algebra matrices T a are to 
be inserted into the Wilson line at the point x, and that we are summing 
over repeated indices. This means that each T a < K\A > is a new Wilson 
line obtained from the original line < K\A > by leaving the form of the 
loop unchanged, but inserting the matrix T a into that loop at the point x. 
In Figure 6 we have illustrated this mode of insertion of Lie algebra into 
the Wilson loop. Here and in further illustrations in this section we use 
Wk{A) to denote the Wilson loop. Note that in the diagrammatic version 
shown in Figure 6 we have let small triangles with legs indicate dx % . The legs 
correspond to indices just as in our work in the last section with Lie algebras 
and chord diagrams. The curvature tensor is indicated as a circle with three 
legs corresponding to the indices of F' a s . 

Notation. In the diagrams in this section we have dropped mention of the 
factor of (l/47r) that occurs in the integral. This convention saves space in 
the figures. In these figures L denotes the Chern-Simons Lagrangian. 

Remark. In thinking about the Wilson line < K\A >= tr(Pe-^ A ), it is 
helpful to recall Euler's formula for the exponential: 

e x = lirrin^ozil + x/n) n . 

The Wilson line is the limit, over partitions of the loop K, of products of 
the matrices (1 + A(x)) where x runs over the partition. Thus we can write 
symbolically, 

< K\A >= J] (1 + A(x)) 
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Figure 6 — Lie algebra and Curvature Tensor insertion into the 

Wilson Loop 

= \{{l+At{x)T a dx k ). 

It is understood that a product of matrices around a closed loop connotes 
the trace of the product. The ordering is forced by the one dimensional 
nature of the loop. Insertion of a given matrix into this product at a point 
on the loop is then a well-defined concept. If T is a given matrix then it is 
understood that T < K\A > denotes the insertion of T into some point of 
the loop. In the case above, it is understood from context in the formula that 
the insertion is to be performed at the point x indicated in the argument of 
the curvature. 

Remark. The previous remark implies the following formula for the varia- 
tion of the Wilson loop with respect to the gauge field: 

5 < K\A > /5(A a k (x)) = dx k T a < K\A > . 

Varying the Wilson loop with respect to the gauge field results in the in- 
sertion of an infinitesimal Lie algebra element into the loop. Figure 7 gives 
a diagrammatic form for this formula. In that Figure we use a capital D 
with up and down legs to denote the derivative 5/5(Al(x)). Insertions in the 
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Wilson line are indicated directly by matrix boxes placed in a representative 
bit of line. 



Figure 7 — Differentiating the Wilson Line 

Proof. 

5 < K\A > /5(A a k (x)) 
= 5U(l + Al(y)T a dy k )/5(Al(x)) 

yeK 

= J] (1 + A a k (y)T a dy k )[T a dx k ] ]J (1 + A a k {y)T a dy k ) 

y<xeK y>xeK 

= dx k T a < K\A > . 

Fact 2. The variation of the Chern-Simons Lagrangian S with respect to 
the gauge potential at a given point in three-space is related to the values of 
the curvature tensor at that point by the following formula: 

F? s (x) = e rst 6S/5(A?(x)). 

Here e a b c is the epsilon symbol for three indices, i.e. it is +1 for positive 
permutations of 123 and —1 for negative permutations of 123 and zero if any 
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two indices are repeated. A diagrammatic for this formula is shown in Figure 



Figure 8 — Variational Formula for Curvature 



With these facts at hand we are prepared to determine how the Witten 
integral behaves under a small deformation of the loop K. 

Theorem. 1. Let Z(K) = Z(S 3 ,K) and let SZ(K) denote the change of 
Z(K) under an infinitesimal change in the loop K. Then 



where Vol = e rst dx r dx s dx t . 

The sum is taken over repeated indices, and the insertion is taken of the 
matrices T a T a at the chosen point x on the loop K that is regarded as the 
center of the deformation. The volume element Vol = e rst dx r dx s dx t is taken 
with regard to the infinitesimal directions of the loop deformation from this 
point on the original loop. 

2. The same formula applies, with a different interpretation, to the case 
where x is a double point of transversal self intersection of a loop K, and the 
deformation consists in shifting one of the crossing segments perpendicularly 
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to the plane of intersection so that the self-intersection point disappears. In 
this case, one T a is inserted into each of the transversal crossing segments so 
that T a T a < K\A > denotes a Wilson loop with a self intersection at x and 
insertions of T a at x+e 1 and x+e 2 where e\ and e 2 denote small displacements 
along the two arcs of K that intersect at x. In this case, the volume form is 
nonzero, with two directions coming from the plane of movement of one arc, 
and the perpendicular direction is the direction of the other arc. 



This completes the formalism of the proof. In the case of part 2., a change 
of interpretation occurs at the point in the argument when the Wilson line 
is differentiated. Differentiating a self intersecting Wilson line at a point of 
self intersection is equivalent to differentiating the corresponding product of 
matrices with respect to a variable that occurs at two points in the product 
(corresponding to the two places where the loop passes through the point). 
One of these derivatives gives rise to a term with volume form equal to zero, 
the other term is the one that is described in part 2. This completes the 
proof of the Theorem. / / 

The formalism of this proof is illustrated in Figure 9. 



Proof. 





(integration by parts and the boundary terms vanish) 




16 




Figure 9 — Varying the Functional Integral by Varying the Line 
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In the case of switching a crossing the key point is to write the crossing 
switch as a composition of first moving a segment to obtain a transversal 
intersection of the diagram with itself, and then to continue the motion to 
complete the switch. One then analyses separately the case where x is a 
double point of transversal self intersection of a loop K, and the deformation 
consists in shifting one of the crossing segments perpendicularly to the plane 
of intersection so that the self-intersection point disappears. In this case, 
one T a is inserted into each of the transversal crossing segments so that 
rparpa < x\A > denotes a Wilson loop with a self intersection at x and 
insertions of T a at x + e\ and x + €2 as in part 2. of the Theorem above. The 
first insertion is in the moving line, due to curvature. The second insertion 
is the consequence of differentiating the self-touching Wilson line. Since this 
line can be regarded as a product, the differentiation occurs twice at the point 
of intersection, and it is the second direction that produces the non-vanishing 
volume form. 

Up to the choice of our conventions for constants, the switching formula 
is, as shown below (See Figure 10). 

Z(K + ) - Z(K_) = (Am/k) J DAe {ik '^ s T a T a < K„\A > 

= (4m/k)Z(T a T a K^), 

where K** denotes the result of replacing the crossing by a self-touching 
crossing. We distinguish this from adding a graphical node at this crossing 
by using the double star notation. 

A key point is to notice that the Lie algebra insertion for this difference 
is exactly what is done (in chord diagrams) to make the weight systems 
for Vassiliev invariants (without the framing compensation). Here we take 
formally the perturbative expansion of the Witten integral to obtain Vassiliev 
invariants as coefficients of the powers of (l/k n ). Thus the formalism of the 
Witten functional integral takes one directly to these weight systems in the 
case of the classical Lie algebras. In this way the functional integral is central 
to the structure of the Vassiliev invariants. 
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Figure 10 — The Difference Formula 



3.1 The Loop Transform 

Suppose that tp(A) is a (complex valued) function defined on gauge fields. 
Then we define formally the loop transform ip(K), a function on embedded 
loops in three dimensional space, by the formula 

$(K) = J DAip(A)W K (A). 

If A is a differential operator defined on ip(A), then we can use this integral 
transform to shift the effect of A to an operator on loops via integration by 
parts: 

A^(K) = J DAAiP(A)W K (A) 

= -J DAiP(A)AW K (A). 

When A is applied to the Wilson loop the result can be an understandable 
geometric or topological operation. In Figures 11, 12 and 13 we illustrate 
this situation with diagrammatically defined operators G and H. 
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We see from Figure 12 that 

Ghp(K) = 8$(K) 

where this variation refers to the effect of varying K by a small loop. As 
we saw in this section, this means that if ip{K) is a topological invariant of 
knots and links, then Gi/j(K) = for all embedded loops K. This condition 
is a transform analogue of the equation Gip(A) = 0. This equation is the 
differential analogue of an invariant of knots and links. It may happen that 
5ip(K) is not strictly zero, as in the case of our framed knot invariants. For 
example with 

ljj(A) = e {ik/4n) / tr(AAdA+(2/3)AAAAA) 

we conclude that Gip{K) is zero for flat deformations (in the sense of this 
section) of the loop K, but can be non-zero in the presence of a twist or curl. 
In this sense the loop transform provides a subtle variation on the strict 
condition Gip(A) = 0. 

In [[TJ] and earlier publications by these authors, the loop transform is 
used to study a reformulation and quantization of Einstein gravity. The dif- 
ferential geometric gravity theory is reformulated in terms of a background 
gauge connection and in the quantization, the Hilbert space consists in func- 
tions ip(A) that are required to satisfy the constraints 

Gtp = 

and 

Hip = 

where H is the operator shown in Figure 13. Thus we see that G(K) can be 
partially zero in the sense of producing a framed knot invariant, and (from 
Figure 13 and the antisymmetry of the epsilon) that H(K) is zero for non- 
self intersecting loops. This means that the loop transforms of G and H 
can be used to investigate a subtle variation of the original scheme for the 
quantization of gravity. This program is being actively pursued by a number 
of researchers. The Vassiliev invariants arising from a topologically invariant 
loop transform should be of significance to this theory. This theme will be 
explored in a subsequent paper. 
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Figure 11 — The Loop Transform and Operators G and H 
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Figure 12 — The Diffeomorphism Constraint 
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4 Wilson Lines, Axial Gauge and the Kont- 
sevich Integrals 

In this section we follow the gauge fixing method used by Frohlich and King 



nj • Their paper was written before the advent of Vassiliev invariants, but 
contains, as we shall see, nearly the whole story about the Kontsevich inte- 
gral. A similar approach to ours can be found in |37|]. In our case we have 



simplified the determination of the inverse operator for this formalism and 
we have given a few more details about the calculation of the correlation 
functions than is customary in physics literature. I hope that this approach 
makes this subject more accessible to mathematicians. A heuristic argu- 
ment of this kind contains a great deal of valuable mathematics. It is clear 
that these matters will eventually be given a fully rigorous treatment. In 
fact, in the present case there is a rigorous treatment, due to Albevario and 
Sen-Gupta [[[[] of the functional integral after the light-cone gauge has been 
imposed. 

Let ( X ■ X - X ) denote a point in three dimensional space. Change to light- 
cone coordinates 



and 



+ 1 i 2 

X = X + X 



x = x 1 — x 2 . 



Let t denote x°. 

Then the gauge connection can be written in the form 

A(x) = A + (x)dx + + A_(x)dx~ + A (x)dt. 



Let CS(A) denote the Chern-Simons integral (over the three dimensional 
sphere) 

CS(A) = (1/4tt) J tr{A AdA+ (2/3)4 AAA A). 

We define axial gauge to be the condition that A_ = 0. We shall now work 
with the functional integral of the previous section under the axial gauge 
restriction. In axial gauge we have that A A A A A = and so 
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CS(A) = (1/4tt) / tr(AAdA). 

Letting d± denote partial differentiation with respect to x ± , we get the fol- 
lowing formula in axial gauge 

A A dA = (A+d-Ao - A d^A + )dx + A dx~ A dt. 

Thus, after integration by parts, we obtain the following formula for the 
Chern-Simons integral: 

CS(A) = (1/2tt) J tr(A + d_A )dx + A dx~ A dt. 
Letting di denote the partial derivative with respect to Xi, we have that 

d + d. = d{ - dl 

If we replace x 2 with ix 2 where i 2 = —1, then <9+<9_ is replaced by 

d\ + d 2 = V 2 . 

We now make this replacement so that the analysis can be expressed over 
the complex numbers. 

Letting 

z = x 1 + ix 2 , 

it is well known that 

V 2 ln{z) = 2ir5(z) 

where S(z) denotes the Dirac delta function and ln(z) is the natural logarithm 
of z. Thus we can write 

= (l/2n)ln(z). 

Note that d+ = d z = d/dz after the replacement of x 2 by ix 2 . As a result we 
have that 

= d+id+d-)- 1 = d + (l/2n)ln(z) = l/2nz. 
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Now that we know the inverse of the operator cL we are in a position to treat 
the Chern-Simons integral as a quadratic form in the pattern 

(-1/2) < A, LA >= -iCS(A) 

where the operator 

L = d_. 

Since we know L^ 1 , we can express the functional integral as a Gaussian 
integral: 

We replace 

Z{K) = J DAe tkCS( > A kr(Pe$K A ) 

by 

Z{K) = J DAe iCS{A kr(Pe$K A /^) 
by sending A to (l/y/k)A. We then replace this version by 
Z(K) = J DAe { - x ^ <A + A> tr{Pe§K A l^ k ). 

In this last formulation we can use our knowledge of L -1 to determine the the 
correlation functions and express Z(K) perturbatively in powers of (l/^/k). 

Proposition. Letting 

< 0(A) >= J DAe^'^^^A)/ J DA e { - l W <A > LA> 

for any functional (f>(A), we find that 

< A a + (z,t)A b + (w,s) >= 0, 

<A a (z,t)A b (w,s) >=0, 
< A a + (z, t)A b {w, s) >= K5 ab 5{t - s)/{z - w) 
where k is a constant. 

Proof Sketch. Let's recall how these correlation functions are obtained. 
The basic formalism for the Gaussian integration is in the pattern 
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< A(z)A(w) >= J DAe(- l IV <A > LA> A(z)A(w)/ J DAe<~- 1 M <A ' LA> 



= ((d/dJ(z))(d/dJ(w))\ J=0 )eW< J > L - 1J > 
Letting G * J(z) = J dwG(z — w)J(w), we have that when 

LG(z) = 5{z) 
(S(z) is a Dirac delta function of z.) then 

LG * J(z) = J dwLG(z - w)J(w) = J dw5(z - w)J(w) = J(z) 

Thus G * J(z) can be identified with L~ x J(z). 
In our case 

G(z) = l/2nz 

and 



L- l J(z) = G * J(z) = J dwJ(w)/(z - w). 



Thus 



< 



J(z),L~ 1 J(z) >=< J(z),G*J(z) >= (1/2tt) J tr(J(z)( J dwJ(w)/(z-w))dz 

= (1/2tt) J J dzdwtr(J(z)J(w))/(z-w). 

The results on the correlation functions then follow directly from differen- 
tiating this expression. Note that the Kronecker delta on Lie algebra in- 
dices is a result of the corresponding Kronecker delta in the trace formula 
tr{T a T}j) = 5 a t > /2 for products of Lie algebra generators. The Kronecker delta 
for the x° = t, s coordinates is a consequence of the evaluation at J equal to 
zero./ / 

We are now prepared to give an explicit form to the perturbative expan- 
sion for 
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,LA> 



DAe (-i/2)<A,LA> tr ( Y[(l + (A/y/k)))/ / DAe^ 1 ' 2 ^ 



LA> 



xeK 



= ^2(l/k n/2 )<f < A(xi)...A(x n ) > . 

„ JK 1 <...<K n 



The latter summation can be rewritten (Wick expansion) into a sum over 
products of pair correlations, and we have already worked out the values of 
these. In the formula above we have written K\ < ... < K n to denote the 
integration over variables xi, ...x n on K so that x\ < ... < x n in the ordering 
induced on the loop K by choosing a basepoint on the loop. After the Wick 
expansion, we get 



Rewriting this in the complexified axial gauge coordinates, the only contri- 



< A(xi)A(aQ > 

=< A%(xi)A%(a/A > T a T a dx + A dt+ < A%{xi)A a + {^ > T a T a dx + A dt 

= (dz-dz')/(z-z')[i/i'\ 




Now we know that 



< A( Xi )A(x'J >=< Alix^A^) > T a T b dx k dx l . 



bution is 



< A a + (z, t)A b (s, w) >= n5 ab 5(t - s)/(z - w). 



Thus 
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where [i/i'} denotes the insertion of the Lie algebra elements T a T a into the 
Wilson loop. 

As a result, for each partition of the loop and choice of pairings P = {xi < 
x'^i = 1, ...m} we get an evaluation D P of the trace of these insertions into 
the loop. This is the value of the corresponding chord diagram in the weight 
systems for Vassiliev invariants. These chord diagram evaluations then figure 
in our formula as shown below: 

m P JK 1 <...<K n i=1 

This is a Wilson loop ordering version of the Kontsevich integral. To see 
the usual form of the integral appear, we change from the time variable 
(parametrization) associated with the loop itself to time variables associated 
with a specific global direction of time in three dimensional space that is 
perpendicular to the complex plane defined by the axial gauge coordinates. 
It is easy to see that this results in one change of sign for each segment of 
the knot diagram supporting a pair correlation where the segment is oriented 
(Wilson loop parameter) downward with respect to the global time direction. 
This results in the rewrite of our formula to 



„ m 

<^>=E(i/nE(-i) lnl ^ / a -<**£)/((* 

m p Jti<...<t n i=1 

where \P j | denotes the number of points (zi,ti) or (z'^ti) in the pairings 
where the knot diagram is oriented downward with respect to global time. 
The integration around the Wilson loop has been replaced by integration in 
the vertical time direction and is so indicated by the replacement of {K 1 < 
... < K n ) with {h < ... < t n } 

The coefficients of l/k m in this expansion are exactly the Kontsevich integrals 
for the weight systems Dp. See Figure 14. 

It was Kontsevich's insight to see (by different means) that these integrals 
could be used to construct Vassiliev invariants from arbitrary weight systems 
satisfying the four-term relations. Here we have seen how these integrals arise 
naturally in the axial gauge fixing of the Witten functional integral. 
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Remark. The careful reader will note that we have not made a discussion of 
the role of the maxima and minima of the space curve of the knot with respect 
to the height direction (t). In fact one has to take these maxima and minima 
very carefully into account and to divide by the corresponding evaluated loop 
pattern (with these maxima and minima) to make the Kontsevich integral 
well-defined and actually invariant under ambient isotopy (with appropriate 
framing correction as well). The corresponding difficulty appears here in 
the fact that because of the gauge choice the Wilson lines are actually only 
defined in the complement of the maxima and minima and one needs to 
analyse a limiting procedure to take care of the inclusion of these points in 
the Wilson line. This points to one of the places where this correspondence 
with the Kontsevich integrals as Feynman integrals for Witten's functional 
integral could stand closer mathematical scrutiny. One purpose of this paper 
has been to outline the correspondences that exist and to put enough light 
on the situation to allow a full story to eventually appear. 
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